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Abstract 

In the present article a family of quantum anharmonic oscillators is studied using 
Hermite's function basis (Fock's basis) in the Hilbert space. The numerical investiga- 
tion of the eigenenergies of that family is presented. The statistical properties of the 
calculated eigenvalues are compared with the theoretical predictions derived from the 
Random Matrix Theory. Conclusions are inferred. 

1 Motivation of the present work 

The simple quantum harmonic oscillator is an excellent model of many physical systems, 
e.g. electromagnetic field. The excitation particles or quanta of the electromagnetic field 
are called photons [1, 2] Moreover it applies for systems of atoms (ions, nuclei) in ideal 
crystals interacting via harmonic attractive force, etc. In the latter case the elementary 
excitation particles of vibrations of crystal lattice or quanta of the sound field are named 
phonons [3, 4]. Both of these quantum fields are bosonic ones [5, 6]. The harmonic potential 
energy operator is only an approximation for the real anharmonic potential energy operator 
of mutual interaction between atoms (ions, nuclei) in real crystals. Therefore the motivation 
of the present paper is an attempt for a more realistic description of quantum anharmonical 
systems. 

2 Quantum harmonic oscillator in D = 3 spatial dimen- 
sions 



Firstly: Quantum harmonic oscillator in D = 3 spatial dimensions is a solvable analytically 
model. In order to make the deliberations easier we reduce our present interest to the 
first quantization case. Therefore the relevant Hilbert space V3 of the complex-valued wave 
functions that are modulus square integrable on the set of {D = 3)-tuples (triples) x 
of real numbers, is isomorphic to a {D = 3)-dimensional tensor (Cartesian) product of the 
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Hilbert spaces Vi of the complex- valued wave functions that are modulus square integrable 
on the set R: ^ 

V3 = L2(R^C,dx) = (g)Vl, (1) 

where 

Vi = L2(R,C,dx,), (2) 

whereas R is the set of the real numbers, and C is the set of the complex numbers. The 
rth harmonic oscillator's eigenfunction ^'r in D = 3 dimensions is, neither symmetrized nor 
antisymmetrized, tensor product of the eigenfunctions in one dimension: 

3 

*r(x) = Y['^rjiXj),:iC^ {xi,...,Xd=3) & R^ T = (n, rD=3) £ N^ (3) 

j=l 

where we used boldface font for the D — 3-tuples x, and r. The Hermite's functions 
(Fock's functions, eigenfunctions of the dimensionless Hamiltonian operator 7i(j) of the quan- 
tum harmonic oscillator in D = 1 spatial dimension) read: 

^rjixj) = Nr^Hr^{xj)exp{-^x^j),Nr. = [V^r,!2'-^]-^ r,- e N, (4) 
where N is a set of natural numbers including zero, whereas: 



Hr,{xj) = (-1)"^ exp(a:,")-^ exjp{-x^A, (5) 



is Vjih Hermite's polynomial [7]. It follows that: 



"1^ 3 3 



*r(x) = NrHri^) exp(--x2) , TV^ = n Nr„Hr(^) - H Hr,(Xj), (6) 



whereas 

h^jf- (7) 



X^ = X • X 



One can also draw a conclusion that the Hilbert space: 

V3 = L2(R^C,dx), (8) 

is composed of the complex-valued wave functions ^ that are modulus square integrable on 
the set R'^. The Hilbert space V3 is separable space, and its orthonormal basis B3 is a set 
of Hermite's functions "i/j. in D — 3 dimensions (Fock's functions in D = 3 dimensions). 
The dimensionless (nondimensional) quantum Hamiltonian operator Ti.^ of the quantum 
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harmonic oscillator in D = 3 dimensions is a sum of quantum Hamiltonian operators 7Y(j) 
of the quantum harmonic oscillators in one dimension: 

^3=1: Tio) = + x'j) = + x^. (9) 

The dimensionless (nondimensional) quantum Hamiltonian operator 7Y(j) of the quantum 
harmonic operator in one spatial dimension is defined as follows: 

T^U)-fj+^l (10) 
The quantum Hamiltonian operator is diagonal in Fock's basis of its eigenfunctions 

Hs^r = er*r, (H) 

and its matrix element ('H3)i,r is equal: 

{n3)l,r = €r5l,r: (12) 

whereas 

3 3 

er = E^-. =E(2^. + 1)' (13) 

is the rth eigenenergy of and where 

Sl,r=l[Sl^,r„ (14) 

is discrete Kronecker's delta in D = 3 dimensions (it is not continuous Dirac's delta in D = 3 
dimensions). The eigenenergies Ej. are simply the sums of all odd natural numbers, and the 
quantum Hamiltonian Ti^ is (direct) sum of diagonal operators in one dimension, and its 
matrix representation is direct sum of diagonal oo x oo matrices (it is poly-index matrix). 

Secondly, let us perform very difficult task consisting of calculating all the flip-flop tran- 
sition amplitudes (hopping amplitudes) from the quantum state xl = x^^r to the quantum 
state *i (s > 0): 

3 poo 

(ms)l,r = (X^)l,r = (*l|i^*r)v3 = 11 / %{x,)x/ ^r,{Xj)dXj . (15) 

The flip-flop transition amplitude (ms)i,r is connected with the processes of emissions and/or 
absorptions of s phonons in D = 3 spatial dimensions, because: 

= Ui^j) - m^^'i^J + ^t)n = [V^(a + a+)]^ (16) 

j=l 3=1 
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where a = (oi, 0,0=3), sl^ — (^1", 0^=3), are the bosonic multiphonon {{D — 3)-phonon) 
annihilation and creation operators in D = 3 spatial dimensions, respectively, and: 



«i*r = V^^'(ri,...,r,-l,...,r^=3)'«/^r = ^Tj + W (ri,...,r,+l,...,ro=3)- (l^) 

It can be easily proven that: 

i^%r=mx;%,r,. (18) 
3=1 

One can calculate the lower transition amplitudes manually, e. g., using recurrence relations, 
matrix algebra, etc., but it is tedious (even for 3 < Sj < 6). The exact formula for all the 
transition amplitudes reads: 

3 

ims)l,r^l[{msj)lj,ry (19) 

i=i 

One can calculate the lower transition amplitudes manually {ms )i , e. g., using recurrence 
relations, matrix algebra, etc., but it is tedious (even for 3 < < 6). If one wants to 
calculate all the transition amplitudes then he must return to the beautiful XIX century 
mathematics methods and after some reasoning he obtains the exact formula [22] : 

i.^^Sj)lj ,rj 

fe/2]h/2] /Ti /ry 

^ ^ ' \^o,^o^^ ^ A,K/,-2A,)! p,!(r,-2p,)! 

.2^+:i-2A,-2,,-i . ^^ 5, + /,+r,-^2A, -2^, + l ^j^ ^20) 

where [•] is entier (step) function, F is Euler's gamma function (compare our result Eq. (20) 
with the formulae in [8, 9]). 

3 Quantum anharmonic oscillator in D = 3 spatial di- 
mensions 

Thirdly, we are ready to investigate the quantum anharmonic oscillator in D = 3 spatial 
dimensions. Its dimensionless Hamiltonian operator 7i|anharm reads: 

^3,a.ha.m = ^3 + E asX^ = 7^3 + ^ ^ ^'^^ [«(.i,..,..=3) " H (%)^^] , (21) 
s=0 (si,...,si3=3)=(0,...,0) j=l 

where S is a D = 3-tuple of degrees of the anharmonicity of the oscillator, and the prefactors 
as are the strengths of anharmonicity. The matrix elements of the anharmonic Hamiltonian 
operator are: 

s s 

(^f,a.iharm)l,r = ^rSl,r + ^ «s(x')l,r = Cr^i^r + ^ ^^^i^s^r, (22) 

s=0 s=0 
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where the representation of the D — 3-dimensional quantum anharmonic oscillator in the 
quantum harmonic oscillator basis is mathematically correct, because the basis is 
a complete set, and the Hilbert space of the eigenfunctions of the anharmonic oscillator 
is isomorphic to the Hilbert space V3 for the harmonic oscillator, provided that the total 
potential energy Uf ^^^^^ of the quantum anharmonic oscillator m. D — Z dimensions: 

W|total(x) = X2 + W|,,haxm(x), (23) 

is bounded from below (there are no scattering eigenstates in D = 3 dimensions) , where the 
anharmonic potential energy W^anharm i^' 

S (Si,- -,513=3) 3 

W|anharm(x) = E «sX^ = E ,..,^0=3) " IK^.)'^]- (24) 

s=0 (si,...,so=3)=(0,...,0) j=l 

It suffices that the D — 3-tuple of degrees of the anharmonicity of the oscillator S = 2S' 
is composed of even numbers and that the strength of anharmonicity ag is strictly positive: 
as > 0, so that Uf^^^^ii;^) oo for |x| — > oo. 

Fourthly, we repeat the "Bohigas conjecture" that the fluctuations of the spectra of the 
quantum systems that correspond to the chaotic systems generally obey the spectra of the 
Gaussian random matrix ensembles. The quantum integrable systems correspond to the clas- 
sical integrable systems in the semiclassical limit [10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21]. 
We emphasize that the "Bohigas conjecture" also holds for the quantum oscillators in D = 3 
dimensions. Having conducted many numerical experiments with different quantum anhar- 
monic oscillators (up to the sextic {D = 3) -dimensional quantum anharmonic oscillators: 
Sj = 6) we draw conclusion that some of them behave like quantum integrable systems, the 
eigenenergies tend to cluster, the histograms of nearest neighbour spacing are closer to the 
Po distribution resulting from the Poisson ensemble, whereas other ones look like quantum 
chaotic systems, their eigenenergies are subject to repulsion, the histograms of NNS are 
closer to the distributions Pi,P2,-Rij derived from the Gaussian Random Matrix ensembles 
[22]. 
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